Memory Effects in the Standard Model for Glasses. 
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The standard model of glasses is an ensemble of two-level systems interacting with a thermal bath. 
The general origin of memory effects in this model is a quasi-stationary but non-equilibrium state 
of a single two-level system, which is realized due to a finite-rate cooling and very slow thermally 
activated relaxation. We show that single particle memory effects, such as negativity of the specific 
heat under reheating, vanish for a sufficiently disordered ensemble. In contrast, a disordered ensem- 
ble displays a collective memory effect [similar to that described by Kovacs for glassy polymers], 
where non-equilibrium features of the ensemble are monitored via a macroscopic observable. An 
experimental realization of the effect can be used to further assess the consistency of the model. 

PACS numbers: 65.60.-|-a, 61.43.Fs 
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Introduction. Low temperature properties of glassy 
and amorphous materials have been an active field of re- 
search for more than 30 years [l[ ; see Q for a review. Ex- 
periments have shown that many characteristics of amor- 
phous materials, e.g., the temperature dependence of the 
specific heat, are universal but different frome those of 
crystals. This evidence has captivated much interest in 
the attempt of producing a coherent theoretical picture 
[H, i]. The two-level system (TLS) model was one of 
the first models to fit the experiments. It soon showed 
to be very successful in describing the low-temperature 
properties of glasses, e.g., the linear temperature depen- 
dence of the specific heat, and gained for itself the def- 
inition of "Standard Model" for glasses [3]. With time 
this model was improved to account for more features of 
amorphous solids [3| and found applications in describing 
low-temperature features of proteins Q . A drawback of 
the model is that there is an excessive freedom in choos- 
ing the distribution of the ensemble parameters. 

Memory effects arise in the model when due to cooling 
down to low temperatures the thermal activation is im- 
peded ^6] . Thus the relaxation time increases to an extent 
that for realistic observation times each TLS is frozen 
in a non-equilibrium, quasi-stationary state, which — in 
contrast to its equilibrium analog — depends on the his- 
tory of the relaxation [6(. Most visible effects of this 
non-equilibrium appear during the subsequent reheat- 
ing, when due to thermal reactivation the single sys- 
tem specific heat becomes negative Q- We shall show 
however that this single particle memory effects do not 
survive the averaging over a sufficiently disordered en- 
semble. In contrast, we propose to implement a memory 
effect, where due to the disorder in the ensemble, locally 
non-equilibrium features of the system are monitored via 
a macroscopic (disorder averaged) observable. This ef- 
fect resembles the one implemented by Kovacs for glassy 
polymers Once the shape of the effect is sensitive 
to dynamic (relaxation) and static (disorder) features of 
the model, its experimental verification would constitute 



a way to further assess the consistency of the model. 
Analogs of the Kovacs effect were recently studied for 
several models of glasses ^ . 




FIG. 1; Upper level occupation n vs. dimensionless tempera- 
ture T/V for a single TLS with e/V = 0.2. Thick curve: n^q. 
The pairs of normal curves refer to cooling (upper curves) and 
reheating (lower curves) with dimensionless rate (from top to 
bottom) a = O.Of , 0.001, 0.0001. 

The standard model of glasses amounts to independent 
particles, each one moving in an asymmetric double- well 
potential with e and V being the energy difference be- 
tween the wells and the barrier height, respectively 
Each particle couples to a thermal bath. The positive 
variables e and V change from one particle to another, so 
that to become observables the single-particle character- 
istics, such as energy or specific heat, should be averaged 
over the joint distribution P{e, V) of e and V. The form 
of P{e, V) is well accounted for in literature 0: 
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where and py are flat distributions with e,-, 

and Vmin <V < Knax- 

There are two regimes in the motion of the single sys- 
tem, i) The thermally activated regime is realized when 
the bath-particle coupling is sufficiently large. At each 
moment of time the particle is then effectively in one 
of the wells, making sudden jumps between them. The 
classical two-state approach is thus a good description of 
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this regime, ii) For low temperatures and weak particle- 
bath couplings only the lowest two energy levels of the 
quantum double well Hamiltonian are relevant and the 
problem reduces to a quantum TLS coupled to a bath 
Here we study only the classical regime. 




FIG. 2: Specific heat C for a single TLS vs. dimensionless 
temperature T/V and e/V = 0.2. Thick curve: equilibrium 
C. Normal curves: C during protocol Q with (from left to 
right) a = 0.0001,0.001,0.01. Dashed curves: the continu- 
ation of each previous cooling protocol by heating with the 
same a (with opposite sign), starting at T/V — 10"^, i.e. 
when n has already relaxed to its zero temperature value. 



Let n and 1 — n be the probabilities for the particle to 
be in the higher and lower well, respectively. Within the 
thermally activated dynamics one has: 
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where e~''*^^+'^^ and e^^^ are the rates of the inter-well 
motion, T — is the bath temperature, and where 70 
is the attempt frequency. Eq. ^ is solved as 
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where r is the relaxation time and ?ie(j(/9) = l/(e^^ + 1) 
the equilibrium value of n reached for i ^ r. At low 
temperatures the relaxation time r becomes very large, 
since there is no enough energy for thermal activation. 
In this regime a freezing temperature T* can be defined 
(see % below which n is essentially frozen-in at T = T* . 
Cooling. Assume that the bath temperature is cooled 
according to the following non-linear protocol: 



(4) 



where a; > is the dimensional cooling rate. This pro- 
tocol is reasonable for low T, since it satisfies the third 
law, not allowing cooling to T = in a finite time. Ex- 
pectedly, for a small rate lo and a high temperature Tf , n 
sticks to its equilibrium value nf.q{j3t) while for lower Tf 
there will not be sufficient time to reach this value, since 
T increases; see Fig. [T] We rewrite ([2|) as: 



a— = [x^ -(- \)n{x) - x^, a = , 

da; 70 



e 



(5) 



where the variable Xt = exp{—f3tV) is introduced and a 
is the dimensionless cooling rate. The solution of ([5]) is: 
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where xq = exp(— y/To). Note from ([5]) that the memory 
about the initial condition Xq is eliminated for xq S> cr. 
If this is satisfied and if a is small, the integral in ([6]) is 
approximated as [a{z) = z^, b{z) = z + z^+^/(/i + 1)]: 
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Eqs. dni E]) leads to the equihbrium value of n: n[t) — 
x^/{l + x^) = neq{/3{t)) This, however, holds under 
neglection of terms a'{x){z — x) and b"{x){z — x)'^ in 
([7|). Thus for the validity of the approximation we 
need: a{x)b'{x) > aa'{x) and b"{x)a < [b'{x)]'^, which 
amounts to aiix^~^ ^ (1 -f x^)"^, and cr/i ^ a;(l -t- x^). 
For /i < 1 and x < 1 we write the relevant conditions as 



a; > (T/x or T > In(cr^)]. 



(8) 



For any finite afi this condition breaks down for low tem- 
peratures a; ^ 0. For these temperatures, x <C fi/i, we 
obtain a non-equilibrium, stationary (time-independent) 
value for n by putting in ([6]) a; = 0. If in addition xq S> cr, 
we put in ((6|) = cxD and get for cr^ ^ 1: 



n(0) 



z^dz 



?I7r+T) = cr^r(l + ^) + o(cr^). (9) 



Compared to n^q, the non-equilibrium n in ([5]) depends 
on the dynamical quantities such as the attempt fre- 
quency 7o and the barrier height V: n is smaller for a 
slower cooling; see Fig. [T] 

Note that the asymmetry fi ^ between the wells 
is crucial for n(0) ^ n^q. For /i ^ we get from the 
integral in Eq. ^ almost equilibrium result: n{0) = 




FIG. 3: Specific heat C versus x/a = a ^^'^ for a single 
TLS with e/V = 0.2. Thick curve: equilibrium C. Normal 
(dashed) curves: C got via cooling (reheating) protocol Q. 
From top to bottom: a = 15 x lO^'', 10"", 0.001. 
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FIG. 4: Average specific heat C versus time t in cool- 
ing from To = WOK (normal curve) to Ti = 0.9K and 
reheating (dashed curve), with two diflterent cooling rates: 
LU = IQ-^iKs)-^ (curves on the left) and lu = 10"^(A's)"^ 
(curves on the right). The probability of the disorder is 
given by © with if < e < 5 iC and 100 K < V < 400 K. 
70 = lO^^s^^, as for experiments by A. Nittke et al. in |2.]. 

1/2 + (^/4) ln[exp(7i;)cr/2] + 0{y?), where 7^ is Euler's 
gamma. 

Specific heat — or the response of the energy en on the 
temperature change — provides more visible effects of the 
memory on the relaxation history. Using ([5]) the equilib- 
rium and the non-equilibrium specific heat are: 

Ceq = edneg/dT = ^?{\nxfx^'[\ + x"]-^ (10) 
C ^en{dt/dT) = ^ia-^x{\nxf [(x^ + l)n{x) - x^"]. 

Since Ceq is zero both for high and low temperatures, it 
displays a maximum at some intermediate temperature; 
see Fig. [2l Under cooling from some high temperatures 
according to ^ , the specific heat Cc shows signs of freez- 
ing: it is smaller than Cgq, saturates quicker to zero, and 
has a smaller maximum. Let us now terminate the cool- 
ing at some temperature T; which is low enough so that 
the energy en relaxed to its stationary value Now 
heat up the bath using the same protocol (|4|) with Tq = T/ 
and a; < 0, and the same dimensionless rate \a\. 

In contrast, the specific heat under heating Ch is seen 
to be negative for sufficiently small temperatures 
This is related to the decrease of the upper-level occupa- 
tion n under reheating; see Fig. [TJ Moreover, Ch ~ — Cc 
at these temperatures: the system keeps memory of the 
cooling stage and still decreases its energy after thermal 
reactivation. Once Ch reaches its negative minimum, it 
quickly increases to the positive maximal value that can 
be larger than the maximum of Ceq- the reheating can 
bring in more thermal instability; see Fig. [5) For higher 
temperatures both Cc and Ch tend to Ceq- The neg- 
ativity of Ch shows that the quasi-stationary state of 
the TLS cannot be viewed as effective equilibrium, as 
far as the reheating is concerned. In order to make the 
meaning of this result more clear, we note that in the 
slow limit, where one decreases \a\ and simultaneously 
increases the time T remains constant, we expect con- 
vergence to equilibrium. Indeed, the temperature region 
where Ch is negative, shrinks to zero as ln(l/|cr|) [see 



([8])], but the magnitude of the negative minimal value 
of Ch in this region does not depend much on a- This 
is seen upon plotting Ch versus x/cr — a^^ exp(— F/T); 
see Fig. [31 However the negativity of Ch , and the very 
difference between Ch and Cc, is sensitive to the values 
e and V of the single-system motion. 

Thus, upon averaging over the disorder — as given 
by ([1]) with experimental values for the parameters of 
P(£, V) — the single-system memory effects gradually 
disappear; see Figs. [H Even though each TLS re- 
mains non-equilibrium, Ch tends to Cc eliminating the 
difference between cooling and heating. The same holds 
for the energy en- 

We shall now discuss another method for displaying 
this non-equilibrium feature. In contrast to the above 
features which are essentially single-system and tend to 
disappear in the presence of disorder, the new method is 
based on the presence of an ensemble. 

Temperature shift protocol. Motivated by Kovacs ex- 
periment we perform the following protocol: 

1. Consider an ensemble of non- interacting TLSs char- 
acterized by a distribution P(e, V)- The ensemble is equi- 
librated at a given high temperature Tq. 

2. Between times t = and te the bath is cooled 
down following The cooling is terminated at a 
low temperature T/ so that the ensemble averaged en- 
ergy {en) = J ededV P{e,V)n{e,V) reached a station- 
ary value. This determines the time tc- Note that (en) is 
observable in experiments using, e.g., heat release mea- 
surements Q- Now {en) equals its equilibrium value: 

{en)U=,^^{e[e''f^ + l]-')^{eneq)- (11) 

This condition defines the temperature T/ = If 
most of the TLSs in the ensemble happen to be described 
at t — tc by a single temperature, then this temperature 
will be close to Tf by definition. Tf turns out to be of 
the same order of the average freezing temperature (T*). 




FIG. 5: The same protocol as in Fig. lU but with stronger 
disorder: OK < £ < 20K and OK < e < 40K (larger curve); 
Lo — 10~^(-R's)~^. In the insert the dependence of the mini- 
mum negative specific heat as a function of the width of the 
energy distribution is displayed for 100 K < V < 400 K and 
lOOK < 1/ < 800 K (dashed line). Emir. = 0. 
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FIG. 6: Ae versus time t for the disordered ensemble with 
probability distribution given by Q . The cooling protocol @ 
started from To = 100 K to Ti = 3K, with lj = 10"^ (Ks)"^ 
Thick curve (stronger disorder): OK < e < lOOK, lOOK < 
V < 800 K, Tf = 16.6 K. Normal curve: OK < £ < 50K, 
lOOK < V < 800 K Tf = 15. 3K. Dashed curve (weaker 
disorder): K < e < 50 K, 100 K < V < 400 K, T/ = 8.8 K. 

3. We want to monitor to what extent the state of 
the ensemble at t — tc is really close to some internal 
equilibrium. To this end, the bath temperature is sud- 
denly switched to Tf, and the resulting evolution of (en) 
is monitored. Due to the sudden switching, the evolution 
is obtained averaging Eq. ^ at the bath temperature Tf, 
and with initial state pT|) : 

Ae = {ent) - {en,q) = (ce"^ (n - l/(e'3^^ + 1))). (12) 

It is seen from (|12p that by our construction Ae should 
be zero both at t = tc and for a very large t — tc- It 
will stay zero for all times t > tc, if the state of (almost) 
each TLS in the ensemble is described by the same tem- 
perature (which need not be equal to that of the bath). 
Yet another case, where Ae is constant for t > tc is 
when there is no disorder in the ensemble. Thus, the 
change of Ae depends both on the disorder and on a 
non-equilibrium state at t = tc- The behavior of Ae for 
experimentally meaningful parameters is shown in Fig. [S] 
Since the change of Ae{t) is finite, a sizable fraction of 
the ensemble is at t = tc far from a local equilibrium. To 
gain more understanding, consider the simplest ensem- 
ble, which is an equal-weight mixture of two TLSs with 
parameters (ei, Vi) and (e2,V2). Eq. (fT2|) impHes 

2Ae = ei(e^^ - e'^)[n{ei,Vi) - l/ie'^f"^ + 1)], 

where n = exp(/3/T/i)/[7o(l -I- exp{-|3fe^))] for i = 1,2 
are the relaxation times (see (O), n(ei,Vi) is given by 
©, and where the temperature /3/ is defined as in pT|) 
summing over the two TLSs ensemble. For the con- 
sidered simplest ensemble, Ae is positive for t > tc- 
This is because the slowest system — e.g., system 1, if 
Ti > T2 — has its non-equilibrium upper-level probability 
7i(ei,Vi) larger than the final equilibrium one [e^f^^ + 
1]^^. In other words, the slowest system is further from 
the equilibrium. The behavior of Ae(i) for an experi- 
mentally relevant disorder distributions ^ is displayed 



in Fig. [6] The fact that e(t) > implies the same expla- 
nation as above: the slow TLSs are further from equilib- 
rium. Two important (and for the present effect general) 
facts seen in Fig. [5] is that the stronger disorder leads to 
i) larger value of Tf and ii) larger maximum of Ae. 

In conclusion, we studied memory effects in the Stan- 
dard Model for glasses. This model, besides describing 
low-temperature properties of many amorphous materi- 
als [2], has important applications in NMR and protein 
physics . It is known from previous works [6,, J] that 
when a single TLS is cooled down to low temperatures, 
the relaxation increases due to impeding of the thermal 
activation, and the system appears in a quasi-stationary, 
non-equilibrium state. In contrast to equilibrium, this 
state depends on the detailed features of the relaxation, 
such as the barrier height or the cooling rate and upon 
reheating manifests itself via a negative specific heat 

We confirmed the latter results by showing that the 
negative magnitude of the reheating specific heat is 
almost insensitive to the decreasing of the cooling- 
reheating rate. Next we showed that the single-particle 
non-equilibrium (memory) effects disappear for a dis- 
ordered ensemble. Since only the latter is experimen- 
tally meaningful, one should question whether the single- 
particle memory effects can be observed at all. Our main 
result is that motivated by Kovacs experiments in 
we designed a protocol which is able to reflect the non- 
equilibrium features of a disordered ensemble. The effect 
is sensitive to the details of the disorder and, if real- 
ized experimentally, it can assess the consistency of the 
model. We have also found two universal features of the 
effect: i) it is more visible for a stronger disorder and ii) 
its sign is determined by the fact that slower elements of 
the ensemble are further from equilibrium. 
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